UNCLASSIFIED 

ad  408  584 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AMD  TECHNICAL  INFORMATION 

CAMERON  STATICN.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used,  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


ACKNOWLEDGMENT 


The  author  sincerely  appreciates  the  helpful  discussions  and  aid 
of  Mr.  A.  Chang  in  the  writing  of  this  note. 


-ii- 


ADDENDUM 


Proof  of  a  General  Relationship 
Used  in  the  Stability  Test  of  Linear  Discrete  Systems 


In  the  stability  criterion  developed  for  the  linear  discrete  systems, 
the  following  general  relationship  is  used  to  reduce  the  stability  constraints. 

Ak=  '  Bk  =  1 1  (Ak+l +  Bk+l)(Ak-rBk-l)  +  (Ak+rBk+l)(Ak-l+Bk-l^  ' 


1 


k  =  2,  3,.  .  .  ,n-l 


(1) 


=  A  A  .  -  B.  -B,  , . 
k-1  k+1  k-1  k+1 


where  A^  +  B^  are  given  in  the  Appendix. 

The  proof  of  Eq.  (1)  is  based  on  the  combined  use  of  the  following 
three  propositions.  The  reader  should  be  familiar  with  Refs.  1  and  2  in 
order  to  follow  the  notations  and  the  details  of  the  proofs. 

Proposition  1: 

3  2k  -  1  <  n,  A^  are  prime  (not  factorable)  polynomials 

in  the  ring  of  polynomials  of  n  variables  a^,  a^,  a^,.  •  •  ,  a^, 

Proof: 

We  readily  notice  that  A^  +  B^  is  prime.  Assume  A^  j  +  B^  j  is 
prime,  we  show  A^  +  is  prime  too.  The  proof  for  A^  -  B^  is  analogous 
to  Ak+Bk- 

Suppose,  on  the  contrary,  that  Aj^+Bj^  is  not  prime.  Let  P  and  Q 
be  polynomials  into  which  ■^jc  +  ®jc  factors, 


AktBk=pQ 


m 


Expanding  the  determinant  (see  Appendix)  for  A^  +  B^  in  cofactors 
of  the  last  row  we  get: 

PQ=  Ak  + Bk=  a0G+ anG*  (3) 

9|C 

where  G  and  G  are  the  cofactors  of  aA  and  a  . 

U  n 
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By  inspection,  it  can  be  seen  that  except  for  a  relabeling  of  some 


indices,  G  with  an=  0  and  G  with  an  =  0  are  of  exactly  the  same  form  as 
+  B, 


Ak-1  '  "k-1 


n  -  0 

(the  restriction  2k  -  2  <  n  is  needed  here  and  for  later  dis- 
cussions  we  require  2k  -  1  <  n).  Therefore  G  I  and  G  I  are 

prime.  'a-=°  '  a"  =  0 

If  we  put  a  =  0  in  (3)  and  divide  by  P, 

a0°|  „ 

Ql 


n 


a  =  0 
n 


(4) 


Also, 


a  G 
n 


a  =  0 
n 


a0=  0 


ao=° 


ao=° 


From  (4)  and  (5)  and  the  fact  that  Q  is  a  polynomial  and  G 
are  prime,  it  follows  that  P  must  be  of  the  form: 

where  P*  is  a  polynomial  and  the  c!s  are  constants. 
Similarly  for  Q  we  can  write 

Q  =  (c  .a  +  c.an  +  c,  +  a  a  Q')< 

4  n  5  0  6  On 


a  =  0 
n 


(5) 


ao=0 


(6) 


(7) 


where  Q'  is  a  polynomial  and  the  c's  are  constants. 

Because  A^  +  B^  is  identically  equal  to  PQ,  all  of  the  constants  in 
(6)  and  (7)  are  necessarily  zero.  This  implies  PQ  a  0  if  a^  or  =  0  which 
is  not  true.  Thus  A^  +  B^  does  not  factor  and  hence  a  prime.  By  similar 
reasoning  A^  -  B^  can  be  shown  to  be  prime  too. 

Proposition  2: 


6.  +a^.  divides  6^+2  i-an-(if2)  ’  1  =  1,2,3 . n-1 


l—  n-i 


(8) 


Proof:  To  prove  Eq.  (8),  first  we  observe  that  the  rules  of  calculating 

6,,  ~  (see  Table  1  of  Ref.  2)  from  the  previous  two  rows  in  the  table  and 

(i+21 

forming  the  sums  +  a'  (i+2)  are  identical  to  those  for  calculating 
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$2  and  —  an-2  from  the  fir8t  two  rows.  Hence,  the  relationship,  as 

far  as  divisibility  is  concerned,  between  an  +  a  and  6,  T  a<2>  i.  exactly 

(i)  ?24.*\  2  —  n”2  7 

the  same  as  that  between  6,  +  a'  \  and  6.  ,,Tr  a"Tr/,f  Therefore,  to 

i—  n-i  i+2  — -  n-(i+2)  #2\ 

prove  (8)  it  is  sufficient  to  show  that  a.  +  *  divides  6,  TV"',  .  This  can 


be  shown  as  follows: 


From  Table  (1)  of  Ref.  (2), 

2  2 

62  =  (a0  ‘  an)  ‘  (aOan-l  ~  *  l*  n)  (9) 

and 

(2)  2  2 

an-2  *  <*0  -  an,(aoV2  ‘  a2an>  *  <a0al "  Vla„><a0a„-1  ‘  ala„>  <1(» 

17) 

In  forming  6-  +  a '  '  ,  we  obtain 
l  —  n 

62‘an-2  =<a0  +  an>  +  l  <"> 

and 

62  +  an-2  *  (a0  -  an»  +2  <12> 

where  and  $  are  polynomials  of  the  variables  a,  ‘s. 

1  .  (2)  * 

Therefore  a»+^  a  divides  6?  a'  and  this  proves  proposition  2. 


Proposition  3: 


Vn»  ^  2k  -  1  <  n  , 


\  +  a<ki 

k  —  n-k 


M  -IN”  V  -IV- 
61  62  63 


6k-2  (Ak-l  £  ^-1* 


=  Ak+1  ±  Bk+1  <13> 


Proof:  To  prove  the  above  relationship  we  introduce  the  following  equation 
which  can  be  obtained  from  Table  1  of  Ref.  2 : 


Vl=  (‘k^n-kXOk  -  <£k  >  ■  2'3 . “-1  <“ 

3 

Furthermore,  Marden  has  proved  the  following  identity  in  terms  of  the 
stability  constants  A^’s  and  B^'s.  * 
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\+r  (Ak+l+  Bk+l)(Ak+l  "  Bk+1)  = 


61  62  63  °k-26k-l 


We  assume  that  (13)  is  valid  for  all  k<  m-1;  by  induction  we  will 
show  that  (13)  holds  for  k  =  m. 

Using  (15)  we  get  for  k  =  m 


(Am+l+  Bm+l)(Am+l  '  Bm+1)  _  ”m- rim-z’/m- ! 

61  62  63 

Using  (14)  and  (15)  for  6  to  obtain 


62  ?5  . 

m-2  m-1 


(6  t  a,m)  )(8  -  a<m>  ) 

/  a  +R  HA  -R  v-  _  m  n-mn  m  n-my _ 

'  m+1  m+1''  m+ 1  m+1'  ..m-1  m-Z-rn^  7Z  T7T  ~  7 

(61  6Z  63  ••6m-Z)(Am-l6l  62  •'6m-3) 


The  above  equation  can  also  be  written  as: 

f  6  +  a(m) 

(A  ,.  +  B  .  )(A  -B  .,)=< - 5 - ^ - 

m+1  m+1  m+1  m+1  <.m-2cm-  3  ,  /A  D  , 

6,  0,  ...o  ,vA_.  i"-o_  i ) 

.  ^  I  c  m-4  m-l  m-l 


6  -a(m> 

_ m  n-m 

6”-28”-3...6m.2(Am_1  +  Bm+1) 


From  this  equation  we  can  readily  establish 
62=A2 

63  =  *  Aj  ,  since  6^  =  Aj) 

*4=  6f  62  A2=  A1  4 
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\+l=  (Ak+l  +  Bk+l)(Ak+l"Bk+l)  =  -k-1  k-2  l.k-3  2  7  (15)+ 

61  62  63  6k-26k-l 

We  assume  that  (13)  is  valid  for  all  k<m-l;  by  induction  we  will 
show  that  (13)  holds  for  k  =  m. 

Using  (15)  we  get  for  k  =m 


(A. 


m+l  +  ^m+l^'V+l  ”  flm+l^  ”  Em-1  cm-2  ,m-3 


■  )(A_ 


-  B 


)  = 


m+1 


Using  (14)  and  (15)  for  6  ..to  obtain 

m-i 


(16) 


6  ,5  . 

m-2  m-1 


(6 


m~Trm^em-3 


(6  ..W](i  -a(m|) 

m  n-m,N  m  n-m ' 


•  5m-2>(Am-l6l  S2  ' 


'6m-3» 


f  1  71 


The  above  equation  can  also  be  written  as: 


(A  XI  +  B  x1)(A  ,.-B  ,.): 

m+1  m+1  m+1  m+1' 


6  +a<m) 

m  n-m 


6?1'26^1'3.  .  .  6m  2(A  -B  .) 

1  L  m-  L  m-i  m-i 


“V 

F, 


J 


6  -  a(m) 

_ m  n-m 


’1 


■V 

F. 


(18) 


+ 


From  this  equation  we  can  readily  establish 


63  =  6^6^  =  Aj  A^  .  since  6j  =  Aj) 

«4=6!  6,  63=A^A2A1A2=Af  A* 
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For  simplicity  we  define  Fj  and  F£  as  indicated  in  the  brackets 

of  (18). 

Next,  we  show  that  Fj  and  are  polynomials: 

(a)  A  ,  +  B  ,  and  [  .  . .  ,  6  _]  are  coprime  (i.  e.  ,  have  no 

m-i  —  m-i  l  2  m-  2 

common  factors). 

By  proposition  1,  Am  ^  +_  Is  prime.  Using  (15)  for  k  =  m-1, 


r  Am“2Am_3  r  1 

l6l  62  -’-tm-zJ 

can  be  written  as  a  product  of  prime  factors  of  the  form  +.  B^»  k<m-l. 
Hence 

[A  ,  +  B  .]  and  [  6?1"2®?1-3 .  .  .  6  ,] 

L  m-1  —  m-lJ  L  1  2  m-2J 

are  coprime. 

(b)  By  assumption,  (13)  is  valid  for  k  <  m-1  and  in  particular,  k=  m-2; 

therefore,  by  noting  (13)  we  readily  ascertain  that  ^  +.  B^  ^  is  a 

divisor  of  , 

c  +  a  (m-2) 

m-2-  n-(m-2  ) 

Now  if  we  apply  proposition  2,  it  is  clear  that  A  .  +  B  divides 
l  \  m-x  —  m—  l 

6mIa'  '  from  (18). 
m  n-m 

2  2 

(c)  Since  A  , ,  is  a  polynomial,  6  -  (a  )  is  divisible  by 

'  m+1  m  '  n-m  7 


r  m-2  m-3  i2 

l6l  62  m-2 


Furthermore, 


A  = 
m 


m 

Pm-2Pm-3 


W1  2  - m-2 

is  a  polynomial,  it  follows  that  a^m^  and  therefore  6  +  a^m^ 

r  7  n-m  rn  —  n-m 

divisible  by 


m  —  n-m 


cm-2cm-3  c 
61  62  ••*6m-2’ 


are 


Items  (b)  and  (c)  in  combination  with  (a)  account  for  all  the  terms  of  the 
denominators  dividing  into  the  numerators  in  F^  and  F^«  Hence,  F^  and 
F^  are  polynomials. 
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From  Proposition  1,  Am+j  +_  ®m+l  arC  Prlme*  Since  FjFj 
in  (18)  is  divisible  by  it  follows  that  Am+^+Bm+^  divides  F^  orF£ 

and  A 


m+1 


Bm+1  also  divides  Fj  or  F^  The  degrees  of  Am+1  +  Bm+1  and 


F^  and  F^  are  equal,  so  A-m+^  +  ®m+l  mu8t  con8tant  multiples  of  Fj  and 
F^-  Therefore,  only  these  two  possibilities  may  exist: 


I. 


II. 


A  + 

B 

=  C,F, 

m+1 

m+1 

1  1 

A 

B 

=  D,  F 

m+1 

m+1 

12 

Am+1  " 

Bm+1 

=  C2F1 

A  + 

B 

=  D  F 

m+1 

m+1 

2  2 

(19) 


(20) 


where  Cj,  D^,  C^,  and  are  constants. 
From  (18)  we  note 


26 


m 


f,(a  rBm  ,)+f,(a  i+Bm  ,)  =  -m'vm  A 

l  m-i  m-1  L  m-1  m-1  r  m-  r  m-  5 


Using  (15),  we  obtain 


(21) 


m-2 


2A  =  F .  (A 
m  I  m- 


l-Bm-l>  +  F2(Am-l+Bm-lli2<A 


m 


■B  )(A  +B  )  (22) 
m  m'  v  ' 


To  show  that  for  and  F^  only  possibility  (I)  exists,  we  substitute 
in  the  above  the  second  form,  i.  e.  ,  Eq.  (20) 

2(A  -B  )(A  +B  )  =  L-  (A  ..-B  ..)(A  ,-B  .) 

m  m  m  m  m+1  m+1  m-1  m-1 

+  —  (A  a1  +  B  , .  )(A  .+B  ,)  (23) 

m+1  m+1'  m-1  m-1'  '  ' 

To  show  the  above  does  not  hold  for  all  m*s,  (1)  we  equate  the  coefficients 
2m 

of  a0  on  both  sides  to  obtain  2  =  1/C_  +  1/D_;  (2)  we  also  equate  the  co- 

efficients  of  a4  on  both  sides.  We  notice  that  for  m-even  ,  the  above 
n  - 

(for  finite  C~  and  D~)  is  not  satisfied  while  possibility  (I)  is  satisfied. 

“  “  2  2(m-l) 

Furthermore,  by  equating  the  coefficients  of  a^a*  ^  "  }  on  both  sides  of 

the  above  equation  in  addition  to  (1)  we  find  that  for  m-odd  the  above  is 


also  not  satisfied  while  using  (19 ) ,  it  is  satisfied.  Therefore  (22) 
is  satisfied  only  by  using  (19)  for  F^  and  F^.  Thus  we  finally 
obtain 


(24) 

We  can  show  that  C^=D^=1  for  possibility  (I)  by  equating  the  coefficients 

of  a  on  both  sides  of  Eq.  (24). 
m 

0  =  +  (4-a  A  ,  -  — —  a  A  ) 

—  n  m-1  n  m-1 

2m 

or  1/C^=1/D^  and  from  equating  the  coefficients  of  a^  on  both  sides  of 
Eq.  (24),  we  get,  2=1/Cj+1/Dj.  Therefore,  Cj=Dj=1.  +  This  completes 
the  proof  of  proposition  3. 

If  we  substitute  for  m  the  variable  k  in  (24), 

Ak=Ak'Bk=  2{^Ak+l+Bk+l)(Ak-rBk-l)+(Ak+rBk+l)(Ak-l+Bk-l)i 


(25) 


Ak-lAk+l'Bk-lBk+l 


which  is  (1)  except  for  the  restriction  on  k  to  be  such  that  2k -1  <  n. 
This  restriction,  which  was  only  used  in  proving  proposition  1  in 
fact,  is  not  necessary.  To  see  that  (25)  holds  for  all  values  of  k  <  n, 
consider  for  the  moment  n  to  be  fixed,  say  n=n^.  We  have  already 
shown  (25)  to  be  valid  for  all  k  such  that  2k-l  <  n^,  and  it  remains 
to  show  its  validity  for  all  other  2k-l  >  n^,  k  <  n^.  For  these  values 
of  k,  consider  (25)  for  n=2nQ.  Since  k  <  nQ,  2k-l  <  n=2nQ  so  (25) 
is  valid  for  these  values  of  k  and  for  n=2nrt.  Put  a~  =a  ,  a_  = 

0  2no  no  2no_1 

=a  ,,  a  ,=a,.  Then  the  resulting  determinants  for  k<n.  are 
n  -1  n^+1  1  ° 

identical  to  those  for  n=n^.  Hence,  (25)  is  valid  for  all  k  <  n^, 
which  proves  the  validity  (1). 

+By  noting  Eqs.  (18)  and  (19),  one  can  establish  that  C.D^=1.  Using 
CjD^=l  in  combination  with  2=1/C^+1/D^,  we  can  readily  establish 
that  Cj=Dj=1. 
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